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ABSTRACT 

Extracting stellar fundamental parameters from SPectro-Interferometric (SPI) data requires 
reliable estimates of observables and with robust uncertainties (visibility, triple product, phase 
closure). A number of fine calibration procedures is necessary throughout the reduction pro- 
cess. Testing departures from centro- symmetry of brightness distributions is a useful comple- 
ment. Developing a set of automatic routines, called SPIDAST (made available to the commu- 
nity) to reduce, calibrate and interpret raw data sets of instantaneous spectro-interferograms at 
the spectral channel level, we complement (and in some respects improve) the ones contained 
in the amdlih Data Reduction Software. Our new software SPIDAST is designed to work in 
an automatic mode, free from subjective choices, while being versatile enough to suit vari- 
ous processing strategies. SPIDAST performs the following automated operations: weighting 
of non-aberrant SPI data (visibility, triple product), fine spectral calibration (sub-pixel level), 
accurate and robust determinations of stellar diameters for calibrator sources (and their un- 
certainties as well), correction for the degradations of the interferometer response in visibility 
and triple product, calculation of the Centro- Symmetry Parameter (CSP) from the calibrated 
triple product, fit of parametric chromatic models on SPI observables, to extract model pa- 
rameters. SPIDAST is currently applied to the scientific study of 18 cool giant and supergiant 
stars, observed with the VLTI/ AMBER facility at medium resolution in the K band. Because 
part of their calibrators have no diameter in the current catalogs, SPIDAST provides new deter- 
minations of the angular diameters of all calibrators. Comparison of SPIDAST final calibrated 
observables with amdlib determinations shows good agreement, under good and poor seeing 
conditions. 

Key words: methods: data analysis - techniques: interferometric - stars: late-type 



1 INTRODUCTION 

The power of optical-infrared interferometry to obtain information 
about the astronomical source morphology (including the angular 
size) is now well-established. To properly determine the source 
properties, the quality of the measurements is an issue, which is 
still the subject of active research. As with any other measuring ap- 
paratus, the absolute calibration of the instrument (including atmo- 
sphere) requires careful attention. Derived from the measurement 
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of the mutual degree of coherence of the incident radiation field, 
on spatial frequencies sampled by the aperture- array configuration, 
the final interferometric observables are non-linear mixes of noisy 
quantities, and of parameter estimates with their own uncertainties. 

In this paper, we propose to revisit and extend the existing 
data processing and calibration methods, in the aim to obtain reli- 
able estimates and robust uncertainties for calibrated measurements 
of visibility and complex triple product. The careful reduction pro- 
cess, described in the present paper, has been elaborated for the sci- 
entific study of a sample of 18 bright cool giant and supergiant stars 
(see Table[T]). Measurements were obtained with the VLTI/AMBER 
facility at medium spectral resolution (^=1500) in the K band, us- 
ing triplets of 1.80-m auxiliary telescopes. Observations have been 
conducted during 15 observing nights between May 2009 and De- 
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Table 1. Science targets of our programme measured with VLTI/AMBER. 
Hipparcos parallaxes are in mas. Last column gives the calibrator(s) associ- 
ated with each science target. 



Name 


Spec. Type 


^Hip 


rriK 


Calibrator(s) 


Of Car 


FOII 


10.6(6) 


-1.3(3) 


77 ColA Eri/HR 3282 


ySCet 


KOIII 


33.9(2) 


-0.3(4) 


77 Get 


QfTrA 


K2II 


8.4(2) 


-1.2(1) 


E TrA/y Lih/o Sgr 


Of Hya 


K3II-III 


18.1(2) 


-1.1(2) 


iHya 


^Ara 


K3III 


6.7(2) 


-0.6(2) 


s TrA/o Sgr 


oi CMa 


K2.5Iab 


0.2(4) 


0.3(3) 


HR2411 


SOph 


M0.5III 


19.1(2) 


-1.2(2) 


7 Lib/6: TrA 


rHyi 


M2III 


15.2(1) 


-1.0(4) 


Of Ret 


o\ Ori 


M3III 


5.0(23) 


-0.7(2) 


HR2411 


crLib 


M3.5III 


11.3(3) 


-1.4(2) 


51Hya 


7 Ret 


M4III 


7.0(1) 


-0.5(3) 


Of Ret 


L2PUP 


M5IIIe 


15.6(10) 


-1.8(1) 


HR 3282/7/ Col 


CETau 


M2Iab-b 


1.8(3) 


-0.9(2) 


40 Ori 


TCet 


M5.5Ib/II 


3.7(5) 


-0.8(3) 


7 SclA Eri 


TXPsc 


C7,2(N0)(Tc) 


3.6(4) 


-0.5(3) 


0Psc 


RScl 


C6,5ea(Np) 


2.1(15) 


-0.1(1) 


L Eri 


WOri 


C5,4(N5) 


2.6(10) 


-0.5(4) 


HR 2113/40 Ori 


TWOph 


C5,5(Nb) 


3.7(12) 


0.5(4) 


Sgr/7 Lib 



cember 2010 (under Belgian VISA Guaranteed Time). The aim of 
the present work is to test the data reduction and calibration soft- 
ware, on data with various qualities, that we started to develop in 
2006. 



2 DEFINING THE BASIC INTERFEROMETRIC 
OBSERVABLES 

The coherent flux Ctj is provided by the measurement of the in- 
stantaneous observables delivered by the AMBER Data Reduction 
Software amdlio^ This quantity traces the sine-like molulated com- 
ponent superimposed on the continuum component, in the observed 
intensity distribution. It is computed using diX^ linear fit of the indi- 
vidual i nterferogranis in th e detector plane, for each spectral chan- 
nel (see lTatulli etaDl200l for details). 

At this point, we define the other interferometric observables: 

(i) the squared flux, which is the product of the photometric 
fluxes Pi and pj, associated with the baseline Bij 



fij = PiPj'^ 



(1) 



(ii) the squared visibility, which is the ratio of the squared mod- 
ulus of the coherent flux to the squared flux 

I |2 



vl = 



■'•J 



(2) 



(iii) the complex bispectrum (IWeigeltll 19771) . which is the prod- 
uct of the 3 coherent fluxes contributing to the baseline triplet 
yBi2,B23,B2ij 

bl23 = C12C23C3U (3) 

(iv) the triple product, which is the ratio of the bispectrum to the 
product of the 3 fluxes contributing to the baseline triplet 

^123 
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Figure 1. Spectral distribution (in K) of the number of non-aberrant data, re- 
maining after box-plot filtering, from an exposure of 300 input frames with 
the calibrator 7 Lib. From top to bottom: baselines AO-DO (32 m), DO-HO 
(64 m), AO-HO (96 m), and triplet AO-DO-HO. Blue lines: good seeing con- 
ditions (seeing angle=0.5 '', coherence time=13 ms). Red lines: poor seeing 
(1.7'', 1.3 ms). 



(v) the closure phase, which is the argument of the complex bis- 
pectrum (or triple product) 



if/i23 = argZ?i23 = arg^i23. 



(5) 



^123 - 



/l 2/23/31 



(4) 



3 REDUCING THE RAW DATA 

To process the temporal series of frames (composing the "expo- 
sures") produced by the AMBER instrument, we use the data re- 
duction software amdlib, which consists of a core library of C func- 
tions plus a high-level interface in the form of a For/ccI plugin (a 
high-level language environment, similar to IDL or Matlab, which 
is in the public domain). 

The standard procedure to extract the basic interferometric ob- 
servables from the AMBER raw data includes a frame- selection 
step, which provides a set of "good" frames according t o a given 
criterion, or a sequence of several selection criteria (Ma lbet et alJ 
2011). The criterion used the most is based on the Signal to Noise 
Ratio (SNR) of the fringe-contrast in each frame. Two selection 
methods are provided in standard by amdlib: 

(i) retaining a given percentage of frames sorted according to 
their SNR values (called the "percentage" method), 

(ii) retaining the frames with SNR values higher than a given 
threshold (called the "threshold" method). 

With both methods, the choice of the optimal value (assumed not 
biasing the final results) of the selection criterion (in percentage 
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Figure 2. Squared visibility at 2.2 yum obtained with amdlih (black solid 
lines + error bars) w.r.t. the percentage of selected frames sorted accord- 
ing to their SNR value, for the calibrator y Lib (baseline AO-DO). Top 
panel: good seeing conditions. Bottom panel: bad seeing. Red lines: val- 
ues obtained with our method (solid line: central value; dashed lines: upper 
and lower limits given by the associated uncertainty). 



of threshold) is done a posteriori, from a sequence of processings 
using different values of the selection criterion (see Millour et al. 
I2OO7 , for details). Two drawbacks of this are: the need for a several- 
steps procedure (to be performed manually), and the risk for biasing 
the final results if the selection criterion is ill-determined. 

In our method, that we have implemented in a specific Yorick 
script added to amdlib, we apply the automatic procedure: 

(i) remove only the aberrant measurements (in squared visibility 
and triple product), /6>r each spectral channel independently. Using 
this specific script results in amounts of useful data larger than the 
amount obtained with the standard procedure, without increasing 
biasing eff'ect (only the lower and upper tails of the input data his- 
tograms are rejected); 

(ii) assign a weight to each individual measurement. After many 
preliminary tests, we decided to choose weights based on the SNR 
of the coherent flux, and on the piston deviation; 

(iii) compute the temporal weighted average of the squared vis- 
ibility and triple product for each spectral channel, with uncertain - 
ties derived from the bootstrap technique invented bv lEfronI (119791) . 

These automatic operations mentioned above, needing neither 
several- step selection nor manual choice of the optimal value of the 
selection criterion, are described in more details in the following 
subsections. 



3.1 Rejecting the aberrant data 

To remove the aberrant measurements, we use a combination of 
physical and statistical basic criteria: 

(i) strictly positive sum and product of the fluxes in the pho- 
© tbd RAS, MNRAS 000.FTlfT4] 
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Figure 3. Histograms of the squared visibilities at 2.2 //m, obtained with 
the calibrator 7 Lib (baseline AO-DO), under good (top panel) and poor 
(bottom panel) seeing conditions. Filled turquoise bars: final histogram after 
aberrant-data rejection (our method). Full black steps: initial histogram (no 
selection). Dashed red steps: final histogram after 50%-frame-selection with 
amdlib. 



tometric channels, since negative flux measurements result from 
shortcomings in the determination of the continuum in each spec- 
trum; 

(ii) coherent-flux SNR higher than unity, since inclusion of data 
with low SNR values reduces the reliability of the estimators of the 
interferometric observables (Mi llour et al.lbOOSi) : 

(iii) squared visibility greater than g 1 - 1 .5 xIQR and lower than 
23 -h 1.5 X IQR ("box-plot filtering"), where Q\ denotes the first 
quartile, Q3 the third quartile, and IQR = (23 - Ql) the interquar- 
tile range (Hoaglin et al. 1983). With triple product data, this third 
criterion applies to the squared visibility for each baseline, and box- 
plot filtering on tan if/ 122 is added to the list of criteria. 

Based on the removal of aberrant measurements at the spec- 
tral channel level, rather than on the selection of the "best" frames, 
our approach provides an amount of useful data larger than the 
amount obtained with the standard frame- selection procedure. For 
data shown in Fig.[T] obtained for a "good" seeing on the calibrator 
target 7 Lib, the number of removed input data varies only slightly 
from one spectral channel to the other (typically less than 3%). Be- 
sides, various amount of data are rejected, according to the seeing 
conditions. 

Using the seeing parameters ^o (seeing angle) and tq (coher- 
ence time), the percentages of rejected data, respective to three 
baselines and two seeing conditions are given in the following lines 
(where "good seeing" refers to sq-O.S '', to=13 ms, and "poor see- 
ing" refers to sq-I.I " , to= 1.3 ms): 

• shortest baseline AO-DO (32 m): 2% (good seeing), to 9% 
(poor seeing); 

• intermediate baseline DO-HO (64 m): 2% to 16%; 

• longest baseline AO-HO (96 m): 2% to 25%; 

• in addition, for the baseline triplet AO-DO-HO, the rejection 
of data amounts 7% and 50%. 

Figure |2] shows the variation of the final squared visibility 
(at l.ljim) produced with amdlib, w.r.t. the percentage of selected 
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Figure 4. Interferometric observables given by our reduction code with the calibrator 7 Lib, before calibration (good seeing). Top left panel: squared visibilities 
for the baselines AO-DO (blue), DO-HO (green), and AO-HO (red) w.r.t. the channel index in K. Top right panel: real (blue) and imaginary (red) parts of the 
triple product for the baseline triplet AO-DO-HO. Because of their low levels, uncertainties in the bottom panels are drawn using histograms, rather than error 
bars (with the same colour code as the top panels). 



frames, based on the SNR and pertaining to the calibrator target 
7 Lib, observed with the shortest baseline AO-DO (32 m) in good 
and poor seeing conditions (top panel and bottom panel respec- 
tively). The squared visibility produced with our method is shown 
for comparison. 

In this example, the position of the plateau of each curve indi- 
cates a value of the selection threshold of 50%, giving the a poste- 
riori determination with amdlib. There we find that our automatic 
method gives a squared-visibility value close to the 50% value, 
without a need for the several- steps procedure, in good and poor 
seeing conditions. Besides, Fig. [3] shows that the standard amdlib 
frame- selection method keeps the extremal V^ values in poor see- 
ing condition, which may bias the result of the calculation of the 
average over the frames contained in each exposure. On the con- 
trary, the box-plot filtering used with our method rejects only the 
lower and upper tails of the V^ histogram, thus leading to estimates 
of the temporal average, more reliable than with the frame- selection 
method. 

Next subsection describes the weighting scheme used to aver- 
age the non-aberrant data over each exposure. Note that selecting 
the frames with the standard amdlib procedure is equivalent to as- 
signing unity weights to the selected frames, and null weights to 
the others. 



3.2 Computing the raw observables 

To reduce the instrumental eff'ects on the data, at sthe pectral chan- 
nel level, we compute in each exposure (300 frames in our exam- 
ple) the weighted average of the squared visibility, and of the com- 
plex triple product. For the visibility, as weight associated with each 
frame we use the ratio of the SNR to the relative excursion of the 
piston, where the relative excursion of the piston itself is the ratio 



of the piston excursion (absolute value) for a given frame, to the 
average over the whole exposure. For the triple product, each indi- 
vidual weight is given by the geometrical mean of the three weights 
associated with the three baselines. The final uncertainties are de- 
rived from the variance of the distribution of the weighted means, 
obtained by random sampling with repl acement of the original se - 
ries of data (direct-bootstrap method) (JEfron & Tibshirani1ll993l) . 
Note that our specific script computes the real and imaginary parts 
of the triple product, from which we extract the calibrated closure 
phase (see Sect. 14.21 ). 

Figure |4] shows an example of squared visibility and triple 
product produced by the reduction script (before calibration) from 
one single exposure of 300 input frames, obtained with the cali- 
brator target 7 Lib, for the three baselines AO-DO (32 m), DO-HO 
(64 m), and AO-HO (96 m). 

In complement to these interferometric observables, we com- 
pute two other quantities used further in the calibration procedure 
(see Sect. 1423]) : 

(i) the weight of each exposure, defined as the product of two 
ratios: one is the average of the SNR to its variance; the other is the 
average of the inverse of the piston excursion (absolute value) to its 
variance; 

(ii) the final number of non- aberrant data, for each spectral 
channel. 



4 CALIBRATING THE DATA 

The calibration process is the key point to obtain accurate estimates 
of the "true" (intrinsic) observables of the science targets. To cor- 
rect the measurements for environmental and instrumental insta- 
bilities, we observe reference stars (the calibrators), with known 
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Figure 5. Spectral distribution of the synthetic model surface flux of y Lib, 
produced by the MARCS + Turbospectrum code in K . Dash-dot line: 
black-body continuous spectrum with reff=4660K. Dash-dot-dot line : En- 
gelke continuous spectrum with the same eff'ective temperature. 
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Figure 6. Spectral distribution of the continuum- subtracted flux of 
7 Lib measured with VLTI/AMBER, after spectral calibration. Dashed 
line: model spectrum (including tellurics). 



angular diameters and independently well-defined brightness dis- 
tributions. 

To calibrate the wavelength-dependent measurements ob- 
tained with VLTI/AMBER, we have developed a library of IDL 
functions, included in the SPIDAST (SPectro-Interferometr ic Data 
Analysis Software Tool) modular software suite (Cr uzalebes et alj 
l2008Ll201Qh . which allows one to: 

(i) link each spectral channel to a wavelength value. We use 
a specific method based on the correlation of calibrator mea- 
sured spectra with synth etic templates given by the MARCS model 
JGustafsson et al.ll2008h : 

(ii) measure and correct for the degradation of the spatial co- 
herence. As calibrators, we use stars with brightness distribution of 
Limb-Darkened (LD) discs, for which angular diameters are given 
by fits of synthetic spectra o n the wide-band spectr ophotometric 
measurements in the infrared dCruzalebes et al.boiob . 



4.1 Computing the spectral shifts 

To reach a high precision in the angular diameter estima- 
tion usin g model fitting, spectral calibration is a critical point 
(see e.g.,lwittkowski et"ani2008l : iDomiciano de Souza et"aDl2008l : 



IStefl et aU201lh . Given the lack of any internal instrumental mod- 
ule for wavelength calibra t ion in the optical setup of AMBER 
(iRobbe-Dubois et all l2007l : IPetrov et al] 120071) . amdlib provides 
calibrated wavelengt h tables, computed from a theoretical polyno- 
mial dispersion law cM erand et al .11201 O r), but with coefficients still 
badly known, leading to wavelength shifts up to ~ 10 pixels in the 
detector plane, at medium spectral resolution (JMalbet et al.ll201lh . 
To correct for this drawback, amdlib shifts the wavelength table us- 
ing the correlation of the measured spectrum with a template table 
containing the telluric lines. 

To improve the precision of the wavelength calibration, ob- 
servers usually compute the coefficients of the dispersion law by 
identifying vis ually some prominent telluric lines in the measured 
spect rum (e.g. jKraus et aLl2009l : IWeigelt et al.l201ll : IOhnaka et al.l 
'2011). We propose an alternative automatic approach, based on 
the cross-correlation of the calibrator measured spectra with syn- 
thetic templates produced by the MARCS + Turbospectrum codes 
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Figure 7. Wavelength correction from the initial polynomial dispersion law 
of AMBER, obtained with SPIDAST (in red) and with amdlib (in blue), for 
each spectral channel of the K band (first observing night). 



(JAlvarez & Plezlll998l : |Plezll2012h . including the tellii r ic line s as 
well. This m ethod was previous ly s uggested bv IPlez (12003 ) for 
Gaia, and by iDecin et al.l (l2004h or iDecin & ErikssonI (l2007h to 
calibrate the SPITZER spectrograms. It is particularly suitable for 
stellar spectra showing easily identifiable spectral features, as it is 
the case with our sample of cool giant calibrator s, showing strong 
CO b ands in the 2. 126-2.474 yum spectral range (IMarti-Vidal et al.l 
l201lh . Figure [5] shows the synthetic spectrum produced by the 
MARCS + Turbospectrum code, for the calibrator y Lib. Some 
reference spectral lines are shown for identification of the corre- 
sponding stellar atmospheric elements. 

For each exposure, obtained with each calibrator, and with 
the working instrumental setting, our automatic spectral calibration 
process contains the following steps: 

(i) apply the heliocentric and systemic radial velocity correc- 
tions, and multiply the synthetic spectrum by the atmospheric and 
instrumental transmittance profiles. Atmospheric transmission data 
for the southern sites ( CTIO, Ch ile) are given by the USAF atmo- 
spheric code PLEXUS JCohen. Wheaton & Megeath 2003); 

(ii) remove the continuum parts of the raw and synthetic spec- 
tra, and normalize the resulting spectra. As a first approximation, 
we estimate the continuum by decreasing the spectral resolution to 
^=40. The main drawback of this method is to produce apparent 
"pseudo-continua" that are low er than the real continuum levels, 
even in almost line-free regions (IRix et al .1120041) . Since our goal is 
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Table 2. Angular diameter values (in mas) of the suitable calibrators, found 
in the CCSL, derived from the final fit of the MARCS-model spectrum, and 
found in the JSDC. The third column gives the eff'ective temperature (in K) 
deduced from the spectral type, used for the model. 



Name 



Spec. Type 



0CCSL 



0final 



0JSDC 



Of Ret 


G8II-III 


4780(230) 


- 


2.54(^) 


2.5(2) 


(P2 Ori 


KOIIIb 


4670(230) 


2.20(2) 


2.263(^) 


2.1(2) 


77 Col 


KOIII 


4660(230) 


2.48(3) 


2.38Q) 


2.7(2) 


AUya 


KOIII 


4660(230) 


- 


2.62(-) 


- 


7 Lib 


KOIII 


4660(230) 


- 


231(]A 


- 


oSgr 


KOIII 


4660(230) 


- 


2.50(-) 


2.4(2) 


L Eri 


K0.5IIIb 


4600(220) 


2.18(2) 


2.35(^,) 


2.5(2) 


6'Psc 


K0.5III 


4580(220) 


2.00(2) 


Ml.) 


2.1(2) 


yScl 


Kim 


4510(220) 


2.13(3) 


2.039(14) 


2.2(2) 


HR2113 


K1.5III 


4440(220) 


- 


2.48(2) 


2.4(2) 


sTrA 


K1.5III 


4440(220) 


2.56(7) 


2.43(2) 


- 


T/Cet 


K1.5III 


4440(220) 


3.44(4) 


3.323(f) 


3.3(2) 


HR 3282 


K2.5II-III 


4330(210) 


2.54(4) 


2.32(7) 




HR2411 


K3III 


4260(210) 


1.90(3) 


i.76(g 


- 


51Hya 


K3III 


4260(210) 


2.28(3) 


2-23(1,) 


2.4(2) 
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Figure 8. Value of the angular diameter of each calibrator of our observing 
sample found in the CCSL, versus the final value deduced fro m SPIDAST, 
fitting the MARCS-model spectra on calibrated templates of ICohen et al.l 
(1999) (red dots), and on IRAS-LRS measurements (black dots). Solid 
lines: ± 10% thresholds. 



only to perform spectral calibration by correlation, this drawback 
has no effect on the final result; 

(iii) divide the spectrum in contiguous sub-windows of the 
same size, and comp ute the wavelength shift by phase-correlation 
JVera & Torresll2008l) in each sub- window. If the spectrum is di- 
vided in n sub- windows, we apply a polynomial law of degree n-l 
to correct for the wavelength shifts. Sub-pixel precision is obtained 
by embedding the cross-power spectrum in the middle of a two- 
times larger array of zeros, before computing the phase-correlation 
function, by inverse Fourier Transform of the fina l cross-power 
spectrum (Guizar-Sicair os. Thurman & Fienup z008); 

(iv) finds the position of the barycenter of the correlation peak 
for each sub-window, computes the coefficients of the interpola- 
tion law (assumed to be polynomial), and corrects the global wave- 
length table; 

(v) calculates the residuals of the wavelength shifts, using once 
again the phase-correlation method on the wavelength-corrected 
spectrum. 

For each observing night, the final wavelength table, associ- 
ated with the working instrumental setting, is obtained by ensemble 
average of the corrected tables (of all exposures with all calibrators 
measured with the same instrumental setting), rejecting the tables 
showing wavelength- shift residuals higher than 8 nm (~5 times the 
nominal spectral resolution). Figure |6] shows the final continuum- 
corrected spectrum, given by the spectral-calibration procedure, 
with the calibrator 7 Lib. The model spectrum drawn with the 
dashed line is produced by the MARCS + Turbospectrum code. 

To show the diflference between the wavelength-calibration 
provided by SPIDAST and by amdlib, we plot in Fig. |7] the cor- 
rections from the initial polynomial dispersion law of AMBER, 
given by the two softwares. While amdlib computes only a unique 
correction value applied over the whole spectrum, our method 
computes a polynomial correction (here parabolic from three sub- 
windows), with sub-pixel deviations with respect to the polynomial 
law, caused by the final ensemble averaging of the corrected wave- 
length tables. The choice of the degree of the polynomial correc- 
tion provided by SPIDAST allows the accuracy of the wavelength 



correction to be adapted to the scientific programme, reaching a 
sub-pixel precision if necessary. 



4.2 Calibrating spectro-interferometric data 

In order to assess and correct for the measurement defects on 
the science targets, proper interferometric calibration (not yet 
fully supported with amdlib) is needed. The Instru ment Response 



luiiy supportea wim amdiw) is neeaea. ine instru ment Kesponse 
Function (IRF), also called "transfer function" (IPerrin bOOSi: 



Ivan Belle & van Belle"2005': 'Boden"2007'; ICruzalebes eTaLll201Qh . 
is derived from observations of calibrators as concomitant as pos- 
sible to the ones of the science targets. 

The standard calibration method is based on the assumption 
that a given "raw" (measured) quantity q of any observable ^, is 
equal to the product of the "true" (intrinsic) quantity Q, multiplied 
by a global degradation factor R^ (in other words the IRF). Such 
a factor cannot practically be modellised with sufficient reliability 
(numerous parameters and some of them remaining unknown). 

To calibrate the quantity related to the science target, a second 
assumption applies, as soon as appropriate conditions (described 
later on) are satisfied and can be then determined by using the rela- 
tion 



qsc 

Qsc 



^cal 



Re 



(6) 



with self explanatory notations, stating that the degradation fac- 
tor is identical for both the calibration source and the science tar- 
get. This assertion is all the more valid that the required conditions 
are fullfilled. The relation here above yields an empirical determi- 
nation of (3sci» via R^, relying on commonplace statistical tools. 
The main condition to satisfy is that the IRF is stable enough, and 
the standard approach to avoid too large variations of the IRF is 
to perform sequences, where observation of calibrator and science 
targets are interlaced and repeated several times, rapidly sampling 
the IRF. So, the actual stability requirement is that the IRF might 
be slowly variable over a given "calibrator- science-calibrator" ob- 
servation sequence. Besides, the calibrator and the science target 
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should be close enough regarding time of observation, angular sep- 
aration, and brightness (in the spectral domain of work), so that 
adjustments of the whole interferometer would not be significantly 
modified from one source to the other. In these conditions, only 
the atmospheric turbulence might cause substantial degradations of 
the measurements. For this latter case, servo-loop systems are to be 
used to reduce turbulence eff'ects. 

Usually, calibrators are found in dedicated catalogs and the 
corresponding sources should meet the desired requirements. How- 
ever, the number of calibrators is necessarily limited and somewhat 
depends on the type of science targets, so that the selected calibra- 
tors might stand out from the ideal conditions: point-like source, 
small angular separation, brightness matching. 

In practice, the limited choice of calibrators makes it necessary 
to accept "partially" resolved targets (van Belle & van Belle 2005), 
angular separations counted in degrees, and brightness mismatch 
amounting some units of magnitudes. For example, a bright cali- 
brator is rarely found in a close angular vicinity of a given bright 
science target. Such a situation nevertheless remains convenient for 
the calibration procedure, since switching between science targets 
and calibrators is fast enough and does not aff'ect the configuration 
of the interferometer, and its response as well. 

This procedure starts with a first selection of calibrators via 
the SearchCal toojj, created by the JMMC working group "Cata- 
logue of calibration sources", providing rough estimates for angular 
diameters, with comparatively large uncertainties. Following this 
fir st selection, the sou rces are searched in the c alibrator c atalogues 
of lBorde et aP (l2002h and lMerand etaP (l2005h . in order to find an- 
gular diameters of better quality. Since some calibrators have not 
a diameter in the published catalogues, it is necessary to make our 
own determinations of those diameters (see Subsect. |4.2.1t . More- 
over, in order to control these determinations and to build a homo- 
geneous set of diameters, we apply our own method to recalculate 
all the diameters for the selected calibrators. The good agreement 
found between our determinations and the ones of the Borde's Cat- 
alogue of Calibrator Stars for Long-Baseline Stellar Interferometry 
(CCSL) (as shown in Tabled and in Fig. [8]), attests to the reliability 
of our specific determinations of angular diameters. 

4. 2. 1 Determining calibrator angular diameters 

To derive the IRF for each spectro-interferometric (SPI) observ- 
able, we firstly need to get a reliable estimate of the angular di- 
ameter of eac h calibrator, an d its associated error. Fitting the stel- 
lar models of lKuruczl (1 1979b on the abs o lutely calibrated spectro- 
photometric templates of ICohen et al] (Il999h , the CCSL con- 
tains 374 calibrators, with limb-darkened angular diameters rang- 
ing from 1 to 3 mas. Since our selected calibrators might have 
no angular-diameter estimate (one third of the calibrators of our 
dataset have no value found in the catalogues), or poorly estimated, 
we have added to SPIDAST various routines for determination of 
the angular dianieter fr om indirect methods, presented in detail in 
ICruzalebes et al.l (120100 . These routines compute the angular diam- 
eter: 

(i) by combining the linear diameter, deduced from absolute lu- 
minosity and eff'ective temperature from the Morgan-Keenan (MK) 
spectral type, with the parallax; 

(ii) by using experimental laws based on the interstellar- 
corrected colour index (surface brightness method); 

^ |www.mariotti.fr/searchcal_page.h-Ei| 
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Figure 9. A-/u map of CLV profiles produced by the MARCS + Turbospec- 
trum code, in the K band, for a KOIII-type star. 
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Figure 10. Theoretical interferometric observables of the calibrator 77 Ceti, 
produced by the LD-MARCS model (in colour), and by the UD model with 
the same value of the angular diameter (in black). Top panel: squared vis- 
ibilities for the baselines Gl-Il (blue), DO-Gl (green), and DO-Il (red). 
Bottom panel: real part of the triple product produced by the models, for 
the baseline triplet DO-Gl -II. For each curve, the thickness gives the un- 
certainty. 



(iii) by scaling synthetic spectra on broadband photometric 
measurements (infrared flux method); 

(iv) by fitting synthetic spectra on infrared spectro-photometric 
measurements (see Sect.|6]). 

Table [2] gives the angular diameter values of our calibrators 
found in the CCSL, as well as our final estimates, derived from the 
fit of MARCS + Turbospectrum synthetic spectra on photometric 
or spectro-photometric measurements. In the last column, we also 
give the estimates found in the JMMC Stellar Diameter Catalogue 
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Figure 11. Temporal variation of the IRF {A=2.2 /um) during one observing 
night. Blue dots: response values averaged over each OB, successively ob- 
tained with the calibrators y Lib, s TrA, and o Sgr. Red dots: response val- 
ues interpolated at the time of observation of the science targets (red dashed 
lines: interpolated response). From top to bottom: response in squared vis- 
ibility for the baselines AO-DO, DO-HO, AO-HO, real and imaginary parts 
of the response in triple product for the triplet AO-DO-HO. 



(JSD C) of Lafrasse et alJ (l20ld). with ac curacies between 7 and 
10% (Delfosse 20041: iBonneau et alJ l2006Y 

Figure [8] shows the comparison of our estimates with the 
CCSL values. The mean difference between our estimates and the 
CCSL is 5%. The mean difference for the three targets with the 
available calibrated templates of Cohen (i.e. 02 Ori, y Scl, and 
T] Cet) is lower: 3.5%. This shows the satisfactory agreement be- 
tween our results and the CCSL, which confirms the reliability of 
our approach for determination of the calibrator angular diameters, 
implemented in SPIDAST. 



4.2.2 Modelling the calibrator visibility and triple product 

Once we determined the angular diameter of each calibrator, and its 
associated uncertainty, we derive the IRF for each SPI observable, 
computing the ratio of the observable measured on the calibrator, 
to the synthetic observable, assumed to represent the "true" calibra- 
tor observable. Each synthetic SPI observable is derived from the 
synthetic coherent flux C^J'"^^^ 

(i) the synthetic flux (Eq. [T]) is given by the synthetic coherent 
flux computed for the null baseline (visibility unity); 

(ii) the squared synthetic visibility (Eq. |2]) by the ratio of the 
squared synthetic coherent flux to the synthetic squared flux; 

(iii) the synthetic bispectrum (Eq. [3]) by the product of the syn- 
thetic coherent flux for the three baselines composing each triplet; 

(iv) and the synthetic triple product (Eq. |4j by the ratio of the 
synthetic bispectrum to the cubed synthetic flux. 



Applying the van-Cittert-Zemike theorem, the model coherent flux 
of each calibrator, assumed to emit a luminous intensity with cir- 
cular symmetry, is derived from the Hankel transform of the source 
intrinsic spectral radiance Lx (in W m~^ yt/m~^ sr~^) 

^2 



-4[ 



Lxir) hl^ri 



b^\ rdr, 



(7) 



where is the stellar angular diameter, r the impact parameter, i.e. 
the distance from the centre of the stellar disc (r = at disc center), 
Bij the length of the projected baseline of the (i,j)-pair of apertures, 
and Jo the Bessel function of order zero. 

Since compa ct photospheres are known to deviate from simple 
uniform discs (see HajianetaL 1998, and references therein), radi- 
ally decreasing from their photometric centre, we use the MARCS 
+ Turbospectrum codes to produce models of Centre-to-Limb Vari- 
ation (CLV) profiles, with input parameters derived from the MK 
spectral type. Figure |9] shows the A-ji map of typical CLV profiles 
given by the mo del, where th e reduced radial coordinate yu is given 
by //= Vl-r2 (■Young|l2003h . 

With "partially" resolved calibrators, the deviation from the 
UD model increases with the baseline length. Figure \T0\ shows 
the squared visibility and the real part of the triple product pro- 
duced by the MARCS and the UD models, for the calibrator rj Ceti 
(0finai=3. 32+0.02 mas), deduced from measurements obtained with 
the DO-Il-Gl baseline triplet (Gl-Il=46m, D0-Gl=69m, and 
D0-Il=79m). In this example, the differences in visibility and 
triple product between the UD model and the LD-MARCS model 
are larger than the uncertainty of the synthetic observable s (derived 
by propagating the uncertainty on the angular diameter), hence the 
justification to use the LD model instead of the UD model. 

4.2.3 Correcting for the degradations of the IRF 

The main source of non- stationary instabilities which affect s the 
fringe formation is the atmospheric-phase turbulence (Roddieg 
1981). This issue is partially solved thanks to the use of fringe- 
tracking servo-loops scanning the fringes with cycle periods 
smaller than the seeing coherence time, defined a s the time over 
which the phase fluctuations remain coherent (Davis & Tango| 
1996; Kellerer & Tokovinin 2007). Thus, the residual instabilities 
are expected to cause only slow drifts of the IRF between obser- 
vations of the science targets and their associated calibrators. In 
order to estimate the response at the time when each science tar- 
get is measured, linear interpolation between the successive mea- 
surements on its calibrator(s) is legitimate (Perrin 2003), provided 
that the "calibrator- science-calibrator" observation sequence, re- 

i )eated several times, rapidly s amples the evaluation of the IRF 
Berger. Dumas & Kaufejl201lh . 

To perform the calibration of the squared visibility and triple 
product data, SPIDAST applies the following procedure, working 
at spectral channel level: 

(i) Compute the measured IRF for each exposure, given by the 
ratio of the calibrator measured observable to the calibrator syn- 
thetic observable (at the same spectral resolution, for the same 
baseline or baseline triplet). Values and uncertainties are derived 
from the 4th-order approximation formulae of Winzer (2000). For 
triple product data, the approximation formulae apply separately to 
the real and imaginary parts of the complex ratio. Since the cal- 
ibrator angular-diameter uncertainties are smaller than 10% (see 
Table [3), the uncertainty of a model observable is simply obtained 
from its first partial derivative w.r.t. the angular diameter, multiplied 
by the angular diameter uncertainty. 
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Figure 12. Histograms of calibrated quantities pertaining to the science target 6 Oph, obtained in good (blue) and poor (red) seeing conditions. From top to 
bottom: squared visibilities with the baselines AO-DO, DO-HO, AO-HO, and closure phase with the triplet AO-DO-HO. Long panels of the first column: dif- 
ference between SPIDAST and amdlib. Small panels of the other columns: associated uncertainties with good (blue: SPIDAST, turquoise: amdlib), and poor 
seeing (red: SPIDAST, pink: amdlib). 



(ii) Average the IRF measurements over each OB. As weight 
associated with each exposure, we use the ratio of the associated 
weight (defined in Sect. 13.21 ) to the variance of the calibrator "raw" 
observable. 

(iii) Determine the IRF at the mean time of each exposure ob- 
tained with the science target, from linear interpolation, or poly- 
nomial fit, on the averaged IRF measurements. Figure [TT] shows 
the temporal variation of the measured and interpolated response in 
squared visibility and triple product, obtained during one observing 
night, with three successive calibrators. 

(iv) Divide the observable measured on the science target for 
each exposure of a given OB, by each interpolated value of the 
IRF obtained for the same OB (using Winzer's formulae). Ap- 
plying a method simi lar to that used with the getCal tool of the 
ICal. Inst. Tech.l (l2008h . a weight is associated with each calibrated 
ratio, which combines information on the observations of the cal- 
ibrator and of the science target (science-calibrator angular sepa- 
ration and observation-time delay, coherence time, number of non- 
aberrant data). 

(v) Compute the weighted average over each OB of the ratios 
used in the calibration, which gives the final calibrated SPI observ- 
able (see Subsect. 14.2.41 ). The argument of the final calibrated triple 
product gives the final calibrated closure phase. 



Although the calibration procedure used by SPIDAST is based on 
the standard calibration method, great care has been taken to pro- 
vide reliable uncertainty estimates, thanks to the use of weights 
tracing the data quality at different levels of the processing. Thus, 
SPIDAST produces large uncertainties on final calibrated observ- 
ables obtained with input data of poor quality. Estimating reliable 
error bars is of crucial importance to get calibrated data, usable for 
the fitting process of parametric models, even on data of uneven 
quality (see Sect.[6l). 



4.2.4 Comparing the final calibrated observables with amdlib 

Although the calibration routine provide d with the amdlib pack- 
age is not yet officially validated (.Malbet et al.ll201ll) . we compute 
the differences in squared visibility and closure phase between SP- 
IDAST and amdlib, and compare them with the associated uncer- 
tainties. 

Figure [12] shows the results obtained on the science target 
6 Oph, with the baseline triplet AO-DO-HO, under good and poor 
seeing conditions. The histograms of the uncertainties show that the 
uncertainties in squared visibility estimated with amdlib are smaller 
than with SPIDAST (whatever the seeing), while it is the contrary 
for the uncertainties in closure phase. 
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Figure 13. Left panels: spectral variation of the imaginary part of the triple product, divided by the spectral average of the real part, for 3 science targets (the 
"spectral mean" quoted above each panel corresponds to the spectral average of the displayed quantity; for the left panels, it thus equals the tangent of the 
global closure phase). Central panels: absolute value of the imaginary part of the triple product, divided by the spectral average of the modulus (for the central 
panel, the "spectral mean" equals the sinus of the CSP). Right panels: absolute value of the real part of the triple product, divided by the spectral average of 
the modulus. 



To explain the disagreement between the uncertainties in 
squared visibility, we remind that SPIDAST takes into account in- 
strumental and environmental conditions in the calculation of the 
uncertainties, while amdlib uses only unweighted variance. Thus, 
we consider as underestimated the uncertainties in squared visibil- 
ity produced by amdlib. 

Regarding the closure phase, it is more difficult to explain 
why the uncertainties from amdlib are larger th an the ones from 
SPIDAST. Indeed, no details are found in M albet et alJ (1201 ih . ex- 
pected to describe the amdlib calculation. SPIDAST computes the 
uncertainties on the real and imaginary parts of the calibrated triple 
product, and derives the uncertainties on the final closure phase 
from them. A deep analysis based on the source code of the cal- 
ibration routine provided by amdlib reveals that amdlib computes 
the uncertainties in the final closure phase, adding the uncertainties 
in closure phase of the science target and of the estimated transfer 
function. Since the closure phase is less stable than the triple prod- 
uct, the final uncertainties in closure phase handled with amdlib are 
larger than with SPIDAST. 

To evaluate the agreement between the two procedures, we 
compare the diff'erences between the final calibrated SPl observ- 
ables they produce, to the associated uncertainties. Since the me- 
dian diff'erences are smaller than the median uncertainties (Fig. [12], 
we conclude that the agreement between the two softwares is satis- 
factory, SPIDAST giving more reliable uncertainties than amdlib. 



5 MEASURING THE DEVIATION FROM CIRCULAR 
SYMMETRY 

The closure phase i/zns couples the phases ip of the Fourier trans- 
form of the source brightness distribution at the three spatial fre- 
quencies /i, /2, and /3 - f\ + fi, probed by the three baselines 
in the following way: j/^123 = ^(/i) + ^{fi) - ^(/s)- In the case of 
a centro- symmetrical brightness distribution on-axis, the phase of 
the Fourier transform is uniformly zero and j/^123 is naturally equal 
to zero. If this source is off'-axis, the phase is linear w.r.t. the spatial 
frequency, so that (fif^) - (fifi) + (fifi), and, subsequently, also 
here if/u3 is equal to zero. Thus, a non-zero phase closure is a clear 
indication, at least qualitatively, for a deviation from circular sym- 
metry. 

In Sect. 15.21 we introduce a new parameter (called the centro- 
symmetry parameter), based on the triple product, more sensitive 
to deviation from centro- symmetry than the closure phase. 



5.1 Global closure phase 

The integration of the real and imaginary parts of the triple prod- 
uct 5Jrue, over the observation spectral band [/tmin;^max], leads 
to the global closu re phase, defined by Ragland et al. ( 2QQQ and 
lTatebeetal](l2006h as 



: atan 






(8) 
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Table 3. CSP values and global closure phases, for the science targets. The 
last column gives the ratio of the associated SNRs. 



Name 


CSP (°) 


^band (°) 


SNRcsp/SNR^ 


(JOph 


0.76(3) 


0.16(4) 


7.0 


Of Car 


0.84(3) 


-0.22(5) 


5.7 


L2PUP 


0.90(4) 


0.39(5) 


3.2 


ySCet 


0.98(4) 


0.32(6) 


4.6 


^Ara 


1.12(4) 


-0.22(6) 


8.1 


aTrA 


1.21(5) 


0.20(7) 


8.3 


aHya 


1.21(9) 


-0.16(11) 


9.4 


TWOph 


2.94(5) 


-2.91(5) 


1.1 


CETau 


3.01(11) 


-0.85(17) 


5.5 


rHyi 


3.25(11) 


2.83(13) 


1.4 


01 CMa 


4.41(197) 


2.88(155) 


1.2 


crLib 


5.12(20) 


-1.90(33) 


4.4 


7 Ret 


5.16(37) 


-4.44(115) 


3.6 


TXPsc 


5.34(25) 


-2.17(32) 


3.4 


oi Ori 


8.10(82) 


7.48(84) 


1.1 


WOri 


10.58(46) 


1.35(79) 


13.3 


RScl 


11.63(59) 


-10.63(67) 


1.2 


TCet 


27.88(51) 


-28.79(41) 
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Figure 14. The CSP versus iI/k, the global closure phase in K, for the sci- 
ence targets. The red inset in the upper right comer enlarges the group of 
points located at the bottom, around the null value of the global closure 
phase, pertaining to the targets 6 Oph, a Car, L2 Pup, /? Cet, ^ Ara, a TrA, 
and a Hya. The dash-dot-dash diagonals trace the relation CSP = \iPkV 



A value of the global closure phase close to zero (less than 1 °, 
for our sample) suggests a high degree of centro- symmetry of the 
brightness distribution, related to the observation spectral band. We 
compute the uncertainty on the global closure phase thanks to the 
direct-bootstrap method, using random sampling with replacement 
of the spectral dataset of the calibrated triple product. The values of 
the global closure phases in K, and their associated uncertainties, 
are given in the third column of Table |3] for each science target. 



5.2 Centro-Symmetry Parameter (CSP) 

When the positive and negative values of 3 5Jrue W along the spec- 
tral domain almost mutually compensate (as seen on Fig. [13] for 
TX Psc and W Ori), the spectral integration produces a nearly null 
global closure phase. As mentioned above, this can be taken as a 
hint for centro- symmetry. However, 3 5Jrue (^) is clearly non-null in 
some parts of the spectrum, and this is a hint for a deviation from 
centro- symmetry. To rule out the contradiction, and so as to con- 
sider this latter possibility, we introduce a new estimator, that we 
call the Centro-Symmetry Parameter (CSP), similar to the global 
closure phase, but using instead the absolute value of 3 5Jrue W in 
the numerator, and the modulus of 5Jrue (^) in the denominator. 



CSP = asin "^^" 






(9) 



As for the global closure phase, a small CSP value (less than 2 °, 
for our sample) suggests a high degree of centro- symmetry. Signif- 
icantly high CSP values, as judged from their uncertainties (boot- 
strap method again), require to use asymmetric brightness distribu- 
tion models, for the fitting on the SPI data. The values of the global 
CSP in the ^-band, and their uncertainties, are given in the second 
column of Table [3] 

To illustrate the diff'erence between the global closure phase 
and the CSP, we plot in Fig. [13] three quantities, for 3 scientific 
targets showing a global closure phase close to zero, 

(i) in the panels on the left is shown the imaginary part of the 



triple product 3 5Jrue, divided by the spectral mearjj of the real part 
(^^rue)- Note that the spectral mean of this ratio is equal to the 
tangent of the global closure phase, as defined in Eq. ([8]); 

(ii) in the panels forming the central column is shown the ab- 
solute value of the imaginary part of the triple product |35Jrue|, 
divided by the spectral mean of the modulus (|=^ruel>- Note that the 
spectral mean of this ratio is equal to the sine of the global CSP, as 
defined in Eq. ([9]); 

(iii) in the panels on the right is shown the absolute value of the 
real part of the triple product |^5Jrue|, divided by the spectral mean 
of the modulus (|5Jruel>- 

The three panels on the left show three typical behaviours for 

^ ^true • 

(i) uniformly close to zero {6 Oph); 

(ii) decreasing symmetrically around zero (TX Psc); 

(iii) increasing symmetrically around zero (W Ori). 

For each of these situations, the integration over the spectrum pro- 
duces a nearly null result (hint for centro- symmetry). However 
the CSP allows to detect deviations from centro- symmetry (see 

Fig. [a. 

If the real and imaginary parts of the triple product are 
wavelength-independent, for the observation spectral range, one 
can show that sin CSP = li/^bandl- If not, no analytical relation link- 
ing the two quantities can be derived from Eqs. ([8j and (|9j, because 
of the integrals in the numerators and the denominators. Figure [14] 
shows the CSP in K displayed versus the global closure phase in the 
same spectral band, for each science target, choosing the OB which 
gives the smallest uncertainty on the CSP. Almost all the stars fall 
along the diagonals of the (j/^band. CSP) diagram, which trace the 
relation CSP = l^^bandl- One (W Ori), however, is flagged as asym- 
metrical with the CSP indicator, but not with </^band- This is precisely 
why we favour the CSP over j/^band. as discussed above. T Cet has 



^ defined as the integral w.r.t. the wavelength, divided by the spectral band- 
width ^X = Amax - /^min 
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Figure 15. Fits on the visibility data of S Oph obtained with the basehnes AO-DO, DO-HO, and AO-HO, for one single OB. Left panel: fit of the UD model; 
right panel: fit of the LD-MARCS model (models in dashed lines). Calibrated data are shown without error bars, for clarity purpose. is the best-fit angular 
outer diameter (in mas),;^^ the reduced chi-squared, and R2 the ajusted coefficient of determination. Bottom panels: residuals (calibrated data - fitted values). 



the largest values for CSP and if/hmd as well, indicating large devi- 
ation from centro- symmetry (may be due to strong asymmetries at 
the wavelengths of the CO bands or even suggestive of the presence 
of a binary companion). The second largest CSP valu e per t ains to 
R Scl, recently revealed as a wide binary by Maer cker et alJ(l2012l) . 
using the ALMA array at millimetric wavelengths. In addition, the 
last column of Table [3] shows that, except for T Cet, the SNR for 
the CSP is comparable or larger than the SNR for the global closure 
phase. 



6 FITTING PARAMETRIC CHROMATIC MODELS 

In order to interpret the final calibrated SPI observables, SPIDAST 
provides a fitting routine, based on the modified gradient -expansion 
algorit hm, very similar to the algorithm invented by 'Le venbera 
(I1944I) , and improved by Marquardt ( 1963). The fit applies on any 
calibrated SPI measurements (to be chosen between visibility, flux, 
coherent flux, bispectrum, triple product, and closure phase), us- 
ing a library of single-component or composite parametric chro- 
matic models, characterised by the Fourier spectrum of their in- 
tensity distribution, and the associated first-order partial derivatives 
w.r.t. the model parameters. Figure [TSlshows examples of visibility- 
fitting results obtained with SPIDAST, using two models (UD and 
LD-MARCS). This fitting routine, applied on spectro-photometric 
data, is used to determine the calibrator angular diameters from the 
MARCS + Turbospectrum synthetic spectra (see Subsect. |4.2.U . 

Note that the LITpro softwar^ developed by the JMMC 
working group "Model fitting", uses a set of elementary 

t ;eometrical and cente r-to-limb darkening functions, as well 
Tallon-Bosc e t al. 2008). However, it does not ofifer, contrary to 
SPIDAST, the possibility to fit stellar-disk models with synthetic 
tabulated radiance data (or exitance, for fits on spectro-photometric 



' |www. jmmc.fr/litpro_page.h£i| 



measurements), such those produced by the MARCS + Turbospec- 
trum code. 

Since the uncertainties of the final calibrated data are not nor- 
mally distributed, the covariance matrix, that comes out of the chi- 
squared fit, cannot be used to infer the parameter uncertainties 
JPress et al.i2007 ). To compute the uncertainties, SPIDAST uses the 
residu al-bootstrap method, described in detail in ICmzalebes et alj 
tOld) : 

(i) "synthetic" datasets are produced from random sampling 
with replacement of the Pearson residuals (diflference between cal- 
ibrated and fitted values divided by the uncertainty of the observed 
value), added to the initial fitted values; 

(ii) fitting the model on these new datasets produces a set of chi- 
squared minima, following a probability distribution, from which 
we extract the boundaries of the confidence interval, with a given 
confidence level; 

(iii) the parameter values associated with these chi-squared 
boundaries give the upper and lower limits of the parameter esti- 
mates, leading to asymmetric uncertainties. 

The whole procedure allows the computing of reliable estimates of 
the parameter uncertainties. 



7 CONCLUSION 

In the present paper, we introduce the new SPIDAST software, de- 
veloped since 2006 with the aim to reduce, calibrate, and interpret 
the visibility and triple product measurements obtained with the 
VLTI/AMBER facility. SPIDAST contains a whole set of modules, 
which can be launched separately or in an automatic batch file, and 
summarised herafter: 

(i) The raw data reduction used by SPIDAST computes the 
weighted average of non-aberrant data, at spectral-channel level, 
providing estimates of the SPI observables using an automated pro- 
cedure, while the method presently used with amdlib selects the 
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"best" frames according to a quality threshold determined a poste- 
riori, after several trials. 

(ii) The wavelength calibration procedure performed by SPI- 
DAST provides spectral shifts following a polynomial law, tracing 
them at the channel level. This is done by computing the cross- 
correlation of the measured spectra of the calibrators with their 
synthetic spectra produced by the MARCS model, while amdlib 
only provides a constant spectral shift over the spectrum, from the 
correlation with the telluric lines. 

(iii) For selected calibrators not included in the calibrator cat- 
alogues, SPIDAST provides several routines for estimation of the 
angular diameter with indirect methods, the most accurate being 
the fit of stellar-atmosphere model spectra given by MARCS on 
spectro-photometric data. The calibrator synthetic observables are 
derived using CLV functions produced by MARCS. 

(iv) To obtain an accurate interferometric calibration (via an au- 
tomated procedure), SPIDAST: (1) divides the calibrator raw data 
with the associated synthetic observables in each spectral channel, 
which gives the instrumental response function in squared visibil- 
ity and triple product; (2) interpolates or fits the response at the 
time of each exposure on the science target; (3) divides the science 
raw data with the interpolated/fitted response, which gives the sci- 
ence calibrated observables for each exposure; (4) computes their 
weighted average over each OB. At each processing step, the un- 
certainties are computed thanks to the bootstrap method applied on 
the weighted means. 

(v) Using the real and imaginary parts of the calibrated triple 
product, SPIDAST measures the deviation from centro- symmetry 
of the brightness distribution of each science target in the obser- 
vation spectral band, thanks to a new parameter, more sensitive to 
asymmetries than the global closure phase. 

(vi) Finally, SPIDAST proposes a complete fitting tool, using a 
set of parametric and chromatic models, and accepting input tables 
of flux/intensity synthetic data. 

Such a careful calibration process of SPI data is a crucial step 
for their trustworthy astr ophysical exploitation, w hich is the topic 
of two associated papers JCmzalebes et al.l2013blld) . Parameter ex- 
traction using non-linear fits of source model, as well as aperture- 
synthesis image reconstruction, need reliable estimates of the cal- 
ibrated observables, with robust uncertainties. Our reduction, cali- 
bration, and fitting routines also apply to any other spectral datasets, 
including spectroscopic data. We made the SPIDAST software pub- 
lico: the source code of any program of our software suite can be 
obtained by sending an e-mail to the first author of the present pa- 
per. 
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